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Abstract 

We consider a system of dynamically-modulated photonic resonator lattice undergoing pho¬ 
tonic transition, and show that in the ultra-strong coupling regime such a lattice can exhibit 
non-trivial topological properties, including topologically non-trivial band gaps, and the associ¬ 
ated topologically-robust one-way edge states. Compared with the same system operating in the 
regime where the rotating wave approximation is valid, operating the system in the ultra-strong 
coupling regime results in one-way edge modes that has a larger bandwidth, and is less susceptible 
to loss. Also, in the ultra-strong coupling regime, the system undergoes a topological insulator-to- 
metal phase transition as one varies the modulation strength. This phase transition has no counter 
part in systems satisfying the rotating wave approximation, and its nature is directly related to 
the non-trivial topology of the quasi-energy space. 

PACS numbers: 42.82.Et, 42.70.Qs, 73.43.-f 
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Creating topological effects [1] for both electrons [2] and photons wm are of significant 
current interests. A powerful mechanism for achieving non-trivial topology in a system is to 
dynamically modulate the system in time. For electrons, such modulation can be achieved 


by coupling with an external e 
Floquet topological insulator 


ec 


romagnetic field, and has been used to create an electronic 


Ibl l2lt. For photons, optical analogue of Floquet topological 


insulators has been demonstrated 


22], Also, time-dependent refractive-index modulations 


can be used to create an effective magnetic field 


23], which can break time-reversal symmetry 


and create a one-way edge state that is topologically protected against arbitrary disorders. 

All previous works on the topological behaviors of dynamically modulated systems have 
considered only the weak-coupling regime where the modulation strength is far less than the 
modulation frequency, and applied the rotating wave approximation (RWA). On the other 
hand, in recent years, the study of light-matter interactions in the ultra-strong coupling 
regime, where the rotating wave approximation is no longer valid, is becoming important 


24 


31 ], It is therefore of fundamental importance to understand topological effects beyond 


the rotating wave approximation. 


In this Letter, we analyze a time-dependent Hamiltonian first proposed in Ref. 23] for 


achieving an effective magnetic field for photons. Unlike Ref. 


23], however, here we focus 


the ultra-strong coupling regime where the rotating wave approximation is no longer valid. 
Experimentally, reaching such ultra-strong: coupling regime is in fact relatively straightfor¬ 
ward with current photonic technology 32]. For this system in the ultra-strong coupling 
regime, we show that the topologically protected one-way photonic edge states can persist 
over a broad parameter range. Compared with the weak-coupling regime, the topologically 
protected one-way edge state is less susceptible to intrinsic losses. We also show that, as 
one varies the modulation strength, there is a topological phase transition that is uniquely 
associated with the ultra-strong coupling regime, and has no counter part in weak-coupling 
systems. 

We start with the Hamiltonian 23] 


H = UJ A ^2 a m a 7n + U B ^ h i h n + ^ V cos (fit + 0 mn )(a^6 n + b ] n a m ), (1) 

m n ( mn) 

which describes a lattice of photonic resonators as shown in Figure [0 The lattice con¬ 
sists of two sub-lattices, each consisting of resonators of resonant frequencies oja and cob, 
respectively, a' (a) and M ( b ) are the creation (annihilation) operators associated with 
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the resonators in the two sub-lattices. The coupling between the resonators are modulated 
dynamically, where V is the maximum coupling strength, hi = ua — is the modulation 
requency, and (p mn is the modulation phase. Such a modulation drives a photonic transition 


33] between nearest neighbor resonators. 



FIG. 1: Lattice composed of two types of resonators A (red) and B (blue). The lines represent 
coupling (or bonds) between nearest neighbors. All coupling strengths are modulated in time 
harmonically. For all bonds along the horizontal direction the modulation phase is zero. The 
bonds along the vertical direction have a spatial distribution of modulation phases as specified in 
the figure. The specified phases are for the hopping matrix elements of the Hamiltonian in Eq. m 
along the positive y direction. 


( 2 ) 


Eq. (JT]) can be transformed to the rotating frame with 

a m (b n ) = a m {b n )e l “ A( ' B '> t . 

The Hamiltonian then becomes 

H = A ~Y^j + 6la m e #m " + alJ) n e i2nt+i(pmn + b ] n a m e~ i2nt ~ i<t,rnn ) , (3) 

(mn ) 

where the first two terms define the Hamiltonian /Lrwa in the rotating wave approximation, 


and the 
Ref. 


as 


; two terms are commonly referred as the counter-rotating terms. 


231 ] considered the weak-coupling regime where V hi and applied the rotating 


wave approximation by ignoring the counter-rotating terms in Eq. ([3]). In this case, H ~ 
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Hrwa becomes time-independent, and is identical to the Hamiltonian of a quantum particle 
on a lattice subject to a magnetic field described by a vector potential A, with <t> mn = 
J/] A ■ dl [34] • Such an effective magnetic field provides a new mechanism for controling 
the propagation of light 351-438]. In particular, it can be used to achieve a topologically 
protected one-way edge state B. 

However, from an experimental point of view, it is important to explore the topological 
behavior of the Hamiltonian of Eq. ([T]) beyond the weak-coupling regime. The required 
modulation in Eq. ([I]) can be achieved electro-optically [32], in which case V ~ 5n/n ■ cuo, 
where n is the refractive index, Sn is the strength of the index modulation, and ujq is the 
operating frequency. For electro-optic modulation of silicon, to minimize free-carrier loss, 
Sn/n ~ 10“ 5 — 10~ 4 [sh]]. Optical communication typically uses an Uq corresponding to 
a free-space wavelength near 1.5 micron. Thus V is typically between 10-100 GHz. On 
the other hand, the modulation frequency O in electro-optic modulation is also on the 
order of 10-100 GHz 32], Thus, experimentally one can readily operate in the regime with 
V ~ fh Similar conclusion can be reached for the acoustic-optical modulation scheme 
implemented in Ref. juj for achieving a photonic gauge potential. Unlike the electronic 


transition, where reaching the ultra-strong coupling regime is a significant challenge 


40 


m, 


for photonic transition [33] it is in fact rather natural that the system operates in the ultra¬ 
strong coupling regime. Thus, systems exhibiting photonic transition can be readily used to 
explore the physics of ultra-strong coupling. 

To explore the topological properties of Eq. ([I]) in the ultra-strong coupling regime, 
we perform a Floquet analysis of the Hamiltonian of Eq. (|3j) . Here we choose a spatial 
distribution of the modulation phase as shown in Figure [I] All bonds along the horizontal 
direction have a zero modulation phase. The bonds along the vertical direction have a 
spatial distribution of modulation phases. In the weak-coupling regime, such a distribution 
corresponds to an effective magnetic flux of 7t/2, or 1/4 of the magnetic flux quanta per 
unit cell. The system is therefore topologically non-trivial in the weak-coupling regime. The 
aim of the Floquet analysis is then to see to what extent such non-trivial topological feature 
persist as one goes beyond the rotating wave approximation. 

Our Floquet band structure analysis follows that of Refs. 44, |45]. The system in Figure 
CD is periodic spatially. Therefore the Hamiltonian can be written in the wavevector space 
(fc-space). For each k point, H in Eq. (J3j) has a period in time T = 1 r/fh Therefore, the 
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solution of the equation id/dt\^>) — H |\&) = 0 in general takes the form |T) = e 


where |$(t + T)) = |$(t)) has a periodicity T in time, and is commonly referred to as the 
Floquet eigenstate, e is the quasi-energy and is defined in the temporal first Brillouin zone 
e G [— 7 r/T, 7 t/T] = [—hi, hi]. The quasi-energies and the Floquet eigenstates can be obtained 


by solving numerically the eigenvalue equation: (H — id/dt)\Q) = £|<f>) 44, [45]]. The quasi 


energy thus obtained as a function of k defines the Floquet band structure. 

As an important subtlety, when performing numerical calculation of the Floquet band 
structures, it is essential to perform a gauge transformation such that the resulting Hamil¬ 
tonian has the smallest possible temporal period. For our case here, the transformation of 
Eq. (|2l) , which is a gauge transformation, serves this purpose. While the Hamiltonians H 
(in Eq. (JT[) ) and H (in Eq. (j3]) ) are equivalent to each other since they are related by the 
gauge transformation of Eq. (J2J) , the temporal periods of H and H are 27 t/ 0 and 7r/fi, 
respectively. A key aspect of topological band structure analysis is to identify band gaps 
that are topologically non-trivial. On the other hand, if one analyze H directly, since the 
corresponding temporal first Brillouin zone is smaller, there is additional band folding along 
the quasi-energy axis, which obscures the band gap. The use of H in Eq. (J3]) is in fact quite 
important for the analysis of the topological aspects of the Floquet band structure. 

We now examine the Floquet band structure of the system. H has a spatial periodicity 
of 4a by 2a along the x and y directions, respectively. Thus, its Floquet band structure has 
8 bands, as we can see in Figure El 

As a comparison, we first consider the band structure of the Hamiltonian //rwa as defined 
by ignoring the counter-rotating terms in Eq. ([3]) . TTrwa has a spatial periodicity of 4a by 
la along the x and y directions, respectively. However, to facilitate the comparison with the 
band structure of H, here we plot the bandstructure of H^wa with the spatial periodicity of 
4a by 2a as well. The resulting RWA band structure, as shown in Figure Ek, thus contains 
8 bands. The bands are two fold degenerate. The bandstructure has the same shape for 
different values of V. (Figure Ek)- 

In the RWA band structure, there are two gaps separating the middle group of four 
bands from the upper and the lower groups, each of two bands, respectively. These gaps are 
topologically non-trivial, as can be checked by calculating the Chern number |20J: 
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(4) 
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FIG. 2: Floquet band structure for the Hamiltonian H of Eq. (J3]) , (a) with RWA from Hrwa (red 
solid line) and V = 0.0217 (blue dashed line), (b) V = 0.517, (c) V = 1.117, and (d) V = 1.517. 


where the summation is over the group of bands, each band indiced by a different a. 

A a = ($„(£, t)|iVg|$ a (M)), (5) 

The Chern numbers for the upper, middle and lower groups of bands are +1, -2, +1, re¬ 
spectively. (As a side note, the middle group of four bands can actually be separated into 
two subgroups each consisting of two bands, separated by Dirac points at e — 0. Each of 
the subgroup has a Chern number -1.) The topological analysis here is consistent with the 
association of an effective magnetic held in this system. The gaps remain open for all non¬ 
zero values of V. Thus with the rotating wave approximation there is no phase transition 
as one varies V. 

Having reviewed the band structure of we now consider the Floquet band structure 

of the full Hamiltonian H. Figure Oi shows the cases of V = 0.0217, the Floquet band 
structure H agrees very well with that of the -Hrwa- 

As one increases V to approximately V > 0.H7, the RWA is no longer adequate to describe 
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the band structure of H. Hence, it is no longer possible to interpret the bandstructure using 
the concept of an effective magnetic held. Nevertheless, the two gaps remain open for V 
ranging from near zero to 1.112 (see Figures [2]). Therefore, in this range of V, the Chern 
numbers for the upper, middle and lower groups of bands must remain unchanged at +1, -2, 
+1, respectively, and hence the topological aspects of the band structure remains the same 
as one goes into the ultra-strong coupling regime. 

As V increases from 0, the bands gradually move away from e — 0, and start to occupies 
more of the temporal first Brillouin zone (See Figure [HJo). At V ~ 12, some of these bands 
reach the edge of the temporal first Brillouin zone. Further increase of V then results in the 
folding back of these bands back into the temporal first Brillouin zone and the closing of 
the gaps as we see in Figured with V = 1.112. No gap is found for larger values of V. We 
see that the increase of V induces a topological phase transition: the system behaves as a 
Floquet topological insulator at small V, and a gapless and topologically trivial “metal” at 
large V. 



F(Q) 

FIG. 3: The projected RWA band structure folded into the temporal first Brillouin zone as a 
function of V. The blue regions correspond to the bands. The white regions are the band gap 
regions. 

The values of V required for achieving such a topological phase transition can be estimated 
by folding the RWA band structure into the temporal first Brillouin zone of [—12, 12] (Figure 
[3]), which predicts that the gap closes at V — 1.112. In comparsion, for the full Hamilontian 
of Eq. ([3]), the gap actually closes at V — 1.1112. Therefore, We see that this topological 
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phase transition is directly related to the non-trivial topology of the quasi-energy space. 
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FIG. 4: (a) The bandwidth of the topologically non-trivial gap as a function of V. (b) The projected 
bandstructure at V = 0.5fh 

A key signature of a topologically non-trivial band gap is the existence of a one-way edge 
state in a strip geometry. We therefore calculate the projected Floquet band structure for 
a strip that is infinite in the ^/-direction and has a width of 21a in the x-direction. The 
projected band structure consists of the quasi energy of all the eigenstates of the system as 
a function of k y . The projected band consists of three groups of bands separated by two 
topologically non-trivial band gaps. We observe the existence of one-way edge mode that 
spans these topologically non-trivial band gaps, as shown in Figure |4 }d with V = 0.5fi. Thus 
the one-way edge state persist in the ultra-strong coupling regime. 

For applications of one-way edge modes to carry information, optimizing the bandwidth 
of such a one-wage edge mode is important. Since the one-way edge mode spans the topolog¬ 
ically non-trivial band gap, the size of such a gap becomes a good measure of the bandwidth 
of the one-way edge mode. In Figure 0k, we plot the size of the topologically non-trivial 
band gap, as a function of modulation strength V. The bandwidth increases with V for 
small V, peaks at V = 0.5fi, and then decreases to zero signifying the topological phase 
transition mentioned above. For a given modulation frequency hi therefore, the bandwidth 
of the one-way edge mode maximizes at the ultra-strong coupling regime. 

The one-way edge mode is topologically robust against disorder-induced back-scattering. 
However, such a mode is still susceptible to intrinsic losses of the materials. For practical 































































FIG. 5: The simulation results with the Hamiltonian (jT]) in a 20x20 resonator lattice (uja = 127r c/a, 
ojb = 0). A point source is placed at the location (1,10). The propagation field profiles are plotted 
at time (a) t = 100 a/c with V = 0.0211 and the source frequency oj s = V; (b) t = 10 a/c with 
V = 0.511 and ui s = V; (c) same as (a) and (d) same as (b) with a loss coefficient 7 = 0.1 c/a at 
steady state. 

application then mitigation of the effect of intrinsic loss is important. For a given modulation 
frequency 11 , operating in the ultra-strong coupling regime results in the one-way edge modes 
that have larger group velocities, and hence are less susceptible to loss, as compared to 
operating in the weak-coupling regime. As an illustration, we consider a finite structure 
described by the Hamiltonian of Eq. ([I]) of a size of 20a by 20a (Figure |5]). We consider two 
systems in the weak-coupling and ultra-strong coupling regime, corresponding to V = 0.0212, 
and V = 0.511, respectively. To probe the systems, we place a point source at the location 
(la, 10 a), and choose the frequencies of the point source to be inside one of the topologically 
non-trivial band gaps. Both systems indeed support one-way edge modes as shown in Figures 
[5ji and b. To study the effect of loss, we further include into Eq. (jl|) an extra term 
—7 where 7 = 0.1 c/a is the loss coefficient. The steady-state field 
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distributions for the same source are shown in Figures (5b and d respectively. We see that the 
photons indeed have a much longer propagation distance in the ultra-strong coupling regime 
with V = 0.5Q (Figure EJ1), as compared to the weak-coupling regime with V = 0.020 
(Figure Eb). 

In summary, we consider a system of dynamically-modulated photonic resonator lattice 
undergoing photonic transition, and show that such a lattice can exhibit non-trivial topo¬ 
logical properties in the ultra-strong coupling regime. From an experimental and practical 
point of view, for the same modulation frequency, operating the system in the ultra-strong 
coupling regime results in one-way edge modes that has a larger bandwidth, and is less 
susceptible to loss, as compared to operating the same system in the weak-coupling regime. 
Our work therefore should provide useful guidance to the experimental quest in seeking to 
demonstrate topological effects related to time-reversal symmetry breaking on-chip 32], We 


also show that in the ultra-strong coupling regime, the system undergoes a topological phase 
transition as one varies the modulation strength. This phase transition has no counter part 
in weak-coupling systems, and its nature is directly related to the non-trivial topology of the 


quasi-energy space. In the context o 


the ultra-strong coupling physics 


24 


recent significant fundamental interests in exploring 


31 


40j-|43|, our work points to the exciting prospect of 


exploring non-trivial topological effects in ultra-strong coupling regime. 
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